Spatial decay of the single-particle density matrix in tight-binding metals: 

analytic results in two dimensions 
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Analytical results for the asymptotic spatial decay of the density matrix p(r,r') in the tight- 
binding model of the two-dimensional metal are presented. In various dimensions D, it is found 
analytically and numerically that the density matrix decays with distance according to the power 
law, p(r, r') oc |r - 



-(D+l)/2 
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The computation of ground-state properties of a 
condensed-matter system from the electronic structure 
depends critically upon the computation of the single- 
particle density matrix (DM). If one possesses the oc- 
cupied eigenstates ipii'"^) of ^ single-particle Hamiltonian 
H, then the density matrix at zero temperature can be 
expressed as p(r,r') = J^t occupied V'^*(r')V'^(I■)• The di- 
agonal element of the DM, p{r, r) , is the charge density. 
The electronic energy may be expressed as Tr(pH), and 
in a simple case the electronic part of the interatomic 
forces may be expressed as Tr(— p9H/9R), where R is 
an atomic coordinate. Viewed in the position representa- 
tion, it is clear that the decay of p{r,r') determines how 
locally one can formulate a calculation of the energy or 
forces. This is of special interest in so-called "quantum 
order-TV" methods in modern first-principles computa- 
tional condensed-matter physics ||. 

The DM provides a means to differentiate between a 
metallic and an insulating state. A considerable body of 
work has been devoted to computing the DM in various 
systems. For insulators, it is well established that, for 
|r - r'l ^ oo, p(r, r') - exp(-7|r - r'|) § |. Recently, 
we have published detailed asymptotic expansions for in- 
sulators in one, two and three dimensions 1^. In metals 
the situation is less clear. Analytic results are available 
for the free-electron gas in any dimensionality {D = 1 — 3) 
and the DM exhibits a power-law decay with Gibbs ring- 
ing (from the abrupt cutoff at the Fermi surface at T = 0) 
1^. Little is known about the "tight-binding" case, ex- 
cept in one dimension where the mathematics is trivial. 
One more "realistic" numerical calculation with a den- 
sity functional Hamiltonian has appeared for Al, which 
produces a DM quite similar to the free-electron gas[0. 
In this Brief Report, we provide analytical asymptotic 
results for the decay of the DM in one dimension and, 
for the first time, in two dimensions for special directions 
on a square lattice. The DM is found to decay gener- 
ally as |r — r'|~(-°"'"^^/^ for large |r — r'|, which coincides 
with the free-electron case. Numerics support this law 
for all three dimensions along various lattice directions. 



We provide detailed expressions for the decay depend- 
ing upon the parameters of the one-band tight-binding 
model and the position of the Fermi level. The analytic 
results are confirmed by direct numerical evaluation of 
the DM. 

Let us consider a tight-binding Hamiltonian defined on 
a lattice, 
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where the orthonormal site basis \i) (one electron or- 
bital per site) spans the Hilbert space of the state vec- 
tors. In the case analyzed below, all the site energies, 
Ei — e, and transfer integrals between nearest neigh- 
bours, tij = t, are constant through the lattice. The 
Bloch functions, |k) = TV-^/^ exp{ik • Rj}|j), with 

dispersion = 2tJ2a cos{kaa) for the simple cubic lat- 
tice (say, with lattice constant a — 1 and e = 0) solve the 
eigenproblem for the Hamiltonian (|I]). 

The object to evaluate is the density-matrix operator, 
p, which can be written in the momentum representation 
as p = |k)/k(k|, with /k being the occupation proba- 
bilities for different eigenstates. For an electronic system 
in thermal equilibrium, these probabilities are the Fermi- 
Dirac factors, /k = /(ek) = (1 exp [(ek - Ai)/r])~\ 
where is the Fermi level and T is the temperature. 
The matrix elements of the density-matrix operator in 
the site basis, pij — {i\p\j), are of special interest for 
obtaining the decay properties of the DM in real space. 
These matrix elements can be written in terms of the ma- 
trix elements of the Green's function operator, Gij{e) = 
(i\{e-nr^\3), as: 
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^ / f{e) J2 - £k) exp{fk • (R, - R,)}de 

(2) 
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f{e)G^J{e + ^0)de 



The problem of the DM decay is thus partly reduced to 
the problem of the evaluation of the off-diagonal elements 
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of the Green's function operator in the site basis, which 
is known to be not at all an easy task (see e.g. Ref. ||] 
and references therein). After some standard algebraic 
manipulations and introducing auxiliary integration (see 
e.g. [^), the expression for Gy = for the simple 
cubic lattice can be recast in the following form: 



G..(e) = -- 
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exp{ize/2t}]^i'^° J,y^(z)dz , (3) 



where the integers Va stand for the Cartesian projections 
of the connection vector Rj — R^, and Jy{z) is the Bessel 
function. 

Starting from this point, we are able to proceed ana- 
lytically further only in the particular case of a square 
lattice {D = 2) along the main diagonal, — Vy = v 
(and in the one-dimensional case as well) , when the inte- 
gral over z can be taken exactly [|lO| : 
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exp{ize /2t}J^{z)dz 



^Q.-i/2(l - 2e2) + ^P._i/2(1 - 2e2 



(4) 



if the energy belongs to the band region, < |e| < 1, 
where e = e/At. The functions P^, and Qi, are the asso- 
ciated Legendre functions of the first and second kind, 
respectively. The expression for the DM can then be re- 
cast as 
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where ep = is the dimensionless Fermi level and the 
zero-temperature case is implied. If the Fermi level lies 
above the band, i.e. ep > 1, all the states are occupied at 
zero temperature and the DM is pi, = (5^0j just reflecting 
the completeness of the basis set. This property, together 
with the even character of the integrand in Eq. ^ , allows 
us to rewrite the expression for the DM (for > 0) in 
the following form: 



Q,_i/2(l-2e2)de 



(6) 



and consider for definiteness only ep > 0. 

The integral in Eq. (||) can be simplified in the asymp- 
totic limit of large j/ — > oo by using the asymptotic ex- 
pression for Qi/-i/2(cos(/)) [l^ . 



Qi/-1/2(C0S(^, _ . 

where ((> — cos^^(l — 2e^), so that 
r0o 
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d0, (8) 
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FIG. 1: The dependence of the absolute value of the DM, 
\pu\, on the lattice index v along the main diagonal [1, 1] in 



the square lattice (Rj 



Rt 



{vxjVy) = (i',u)) for a tight- 



binding model of a crystalline metal at zero temperature and 
various positions of the Fermi-level, ep = /i/4t, as marked in 
(a)-(c). The open circles represent the exact numerical result 
obtained from Eq. (^. The crosses and pluses correspond 
to the approximate results obtained by using Eqs. ^ and 
(PH, respectively. The straight solid lines show the power- 
law dependence, p„ oc u~^^^ . 



with 00 ~ cos^^(l — 2e|). The final step revealing the 
explicit asymptotic dependence of the DM on i/ can be 
made if the Fermi level lies not far from the midband 
point and (j>Q <^ 1. In that case, the first term of the 
integrand in Eq. (0) can be expanded in (j), resulting in 
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with '00 = \/i^0o- AH the non-trivial dependence on ly is 
now in the upper limit of the integral, the latter being 
proportional to the Frcsncl integrals C{x) and S{x), 
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and in the asymptotic limit, yj2v(^QpK — s- oo |10 
1 
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As follows from Eq. (|l l|) , the DM decays according to 
the power law, pij oc |Rj — R^l^^^^, at least along the 
main diagonal in the square lattice. 

All the analytical results presented above can be ver- 
ified by direct numerical analysis. In Figs, ^a)-(c), we 
show the dependence of the DM versus v along the diag- 
onal [1,1] in the square lattice for different positions of 
the Fermi level. The exact numerical results (open cir- 
cles) have been obtained by both direct summation over 
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FIG. 2: The dependence of the absolute value of the DM, 
|pi/|, on the lattice index i/: (a) along the direction [1, 0] in the 
square lattice (Rj — Rj = {fxjO) = (i^, 0)) and (b)-(c) along 
the direction [2, 1] in the square lattice (R^— R-/ = {^Vy, Vy) = 
{2v,v)) for the same model as in Fig. gj. The results are 
obtained numerically from Eq. (jij). The straight solid and 
dashed lines show the power-law dependences, oc 
and pi, oc , respectively. 



tlie first Brillouin zone in Eq. (|^) (over 10* points) and 
by integration of Eq. (^) (both methods give identical 
results). The approximate analytic results according to 
Eqs. (^) and ( pT| ) are given by the crosses and pluses, 
respectively. Good agreement between the exact and ap- 
proximate dependencies is evident, even for the relatively 
large values of the Fermi-level position far away from the 
midband region (see Fig. |l|(c)). The solid straight lines 
in Fig. |l| corresponding to the power law, \pv\ oc v~" 
confirm the same law for the DM decay. 
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In order to verify the power-law decay of the DM along 
other directions in the simple square lattice, we have cal- 
culated the DM numerically for these directions and have 
found the same asymptotic behavior, p^, oc v~^l'^ . The 
results for directions [1,0] and [2,1] and different posi- 
tions of the Fermi level are presented in Figs, ^(a)-(c). 

For the special case that the Fermi level lies exactly at 
the band center, ep = 0, note that for the main diagonal 
[1,1], the DM vanishes for > (see Eq. (|)). This 
behavior is also manifested in other directions: for ep ~ 
from pn oc v~^l'^ to oc (see the dashed lines in 
Figs. 11(a)- (b)). Fig. ||(c) shows how this occurs for the 
direction [2, 1] when ep 0. 

The analysis of the DM decay in different dimensions 
can be performed analytically for D = 1 and numerically 
for Z) = 3. First, we look at the one-dimensional system. 
In the case of zero temperature, the integrals in Eq. (|^) 
can be taken exactly, resulting in the following expression 
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FIG. 3: The dependence of the absolute value of the DM, 
jp^j, on the lattice index v. (a) along the direction [1, 0, 0] in 
the simple cubic lattice (Rj — Rj = {^^,0,0) = (z/, 0, 0)) and 
(b)-(c) along the direction [1, 1, 1] in the simple cubic lattice 
(Rj— Rj = (i/^ci'yjVz) = for a tight-binding model of 

a crystalline metal at zero temperature and various positions 
of the Fermi level, ep — p/8t, as marked. The results are 
obtained numerically from Eq. The straight solid and 

dashed lines show the power-law dependences, pi, oc and 
Pv oc respectively. 



for the DM (|^/2t| < 1): 
1 



Pi, — — sm 



/ . -1 A* , 

A"^ 2^+2 



(12) 



This expression has the correct limits for a fully occupied 
(p = 2t) and an empty (p = —2t) band, viz. p^ = 6^q 
and pi, ~ Q, respectively, and shows the power-law decay, 
Pi, oc v^^ . Bearing in mind the decay law, p^, oc 
found for D — 2, we can infer that the generalized law in 
all dimensions is: 



Pi] ^ Pi, cciy 



oc |R, - RJ 
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(13) 



In order to check Eq. (|13|) for Z? = 3, we have cal- 
culated numerically the DM for the simple cubic lat- 
tice along different directions for various positions of the 
Fermi level. The results are shown in Figs, ^(a)-(c), from 
which it is clear that indeed the DM asymptotically sat- 
isfies Eq. ( p^ ) and pi, oc 

Similar to the 2D-case, the mid-band location of the 
Fermi energy brings additional symmetry to the problem 
which can change the asymptotic behaviour of the DM 
along certain directions. For example, if ep = then 
Pi, = for v > along the direction [1, 1,0] and pi, oc 
v^^^"^ along the [1,1,1] direction (see the dashed line in 
Fig. ||(b); Fig. ||(c) demonstrates how the new asymptotic 
behavior appears when the Fermi level approaches the 
midband position). 



4 




FIG. 4: Contour plot of the real-space density matrix for Al 
calculated in the {100} plane for the conventional cubic unit 
cell (the X — y axes are parallel to the bonds). 




FIG. 5; Spatial decay of the real-space density matrix for 
Al calculated for different angular directions, as marked (zero 
angle refers to the bond direction), in the {100} plane for 
the conventional cubic unit cell. The solid circles represent 
calculated values of the density matrix for the free-electron 
gas model with the same electron density as for Al. The solid 
line shows an dependence. 



The above analytical and numerical results have been 
obtained for very simple tight-binding models. This of 
course leaves open the question about the generality of 
our findings. In order to answer this question, at least to 
some extent, we have calculated the single-electron DM 
for a realistic model of fee aluminum (500 atom super- 
cell with the box side of 20.25rA) using an approximate 
density functional Hamiltonian in the local density ap- 
proximation (see Ref. for more detail). The results 
are presented in Figs, ^and ^ The real-space contour 
plot for the DM in the {100} plane for the conventional 
cubic unit cell is shown in Fig. |^. From this plot, we can 
see the isotropic metallic nature of the bonding, in con- 
trast to the case of semiconductors with covalent bonding 
(cf. Fig. 3 in Ref. 0). The spatial decay of the DM along 
different directions in the same symmetry plane for Al is 
shown in Fig. ||, together with the data calculated for 
the free-electron gas model [|j (with the electron density 
being the same as that for Al, i.e. 0.185rA~^). It is 
clearly seen that the DM for Al decays in a very sim- 
ilar fashion to that for the free-electron gas model, i.e. 
p cx r~^. Therefore, these results support the generality 
of our model calculations. 

In conclusion, we have presented analytical and nu- 
merical arguments supporting the power-law decay of the 
density matrix, pij cx |Rj — Ri|~(-°+^^/^, in tight-binding 
models of metals in different dimensions at zero temper- 
ature. The main result is the analytical asymptotic de- 
pendence of the density matrix versus distance along the 
main diagonal in the square lattice (see Eq. (|ll|)). Ap- 
parently, the sharp cutoff induced by the Fermi-Dirac 
distribution at zero temperature in the integration over 
the energy spectrum, independently of the shape of the 
density of states, results in the power-law decay of the 
density matrix in crystalline metals. 
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